Scaling theory of transport in complex networks 
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Transport is an important function in many network systems and understanding its behavior on 
biological, social, and technological networks is crucial for a wide range of applications. However, 
it is a property that is not well-understood in these systems and this is probably due to the lack of 
a general theoretical framework. Here, based on the finding that renormalization can be applied to 
bio-networks, we develop a scaling theory of transport in self-similar networks. We demonstrate the 
networks invariance under length scale renormalization and we show that the problem of transport 
can be characterized in terms of a set of critical exponents. The scaling theory allows us to determine 
the influence of the modular structure on transport. We also generalize our theory by presenting 
and verifying scaling arguments for the dependence of transport on microscopic features, such as 
the degree of the nodes and the distance between them. Using transport concepts such as diffusion 
and resistance we exploit this invariance and we are able to explain, based on the topology of the 
network, recent experimental results on the broad flow distribution in metabolic networks. 



I. INTRODUCTION 

Transport in complex networks is a problem of much 
interest in many aspects of biology, sociology and other 
disciplines. For example, the study of metabolic fluxes 
in organisms is crucial for a deeper understanding of 
how the cell carries its metabolic cycle The use of 
metabolic flux analysis can provide important cellular 
physiological characteristics using the network stoichiom- 
etry and predict optimal flux distributions that satisfy a 
defined metabolic objective. Similarly, information flow 
between the molecules of a biological network provides 
insight for both the network structure and the functions 
performed by the network. Such an example is the con- 
cept of the 'diffusion distance' in a protein-protein inter- 
action network which is used to predict possible inter- 
actions between proteins, simply by studying diffusion 
in the existing network In food webs, energy trans- 
fer between different levels of the web is crucial for the 
organism survival, while spreading of a disease between 
different organisms may affect the regular operation of 
the equilibrated system. Moreover, applications of trans- 
port in complex networks extend to a plethora of other 
systems, such as, for instance, transport of information 
in the Internet, spreading of diseases and/or rumours in 
social networks, etc. Despite its significance, the laws 
of transport in such a complex substrate are yet unclear 
compared to transport in random media [3, . This is 
due to the complexity added by the heterogeneous degree 
distribution in such networks. 

We study transport in real-world biological networks 
and via a model, which possess both self-similar proper- 
ties and the scale-free character in their degree distribu- 
tion. We explain our results with theoretical arguments 
and simulation analysis. We use approaches from renor- 
malization theory in statistical physics that enable us to 
exploit the self-similar characteristics of the fractal net- 
works and develop a scaling theory of transport, which 
we use to address the effects of the modularity and the 



degree inhomogeneity of the substrate. 

Due to the existence of a broad degree distribution, 
transport on a network is different when it is between 
two hubs with a large number of connections k or be- 
tween low-degree nodes. We therefore characterize the 
transport coefficients by their explicit dependence on fci , 
/c2 , and where fci and k2 denote the degree of two nodes 
{ki > fc2), separated by a distance £ (distance is measured 
by the minimum number of links, i.e. it is the chemical 
distance). We study the diffusion time T{£;ki,k2) and 
the resistance R{£', ki, k2) between any two nodes in the 
system. The dependence on ki and k2 is not significant 
in homogeneous systems, but is important in networks 
where the node degree spans a wide range of values, such 
as in biological complex networks. In fact, this depen- 
dence is critical for many other properties as has been 
already shown for e.g. fractality, where traditional meth- 
ods of measuring the fractal dimension may fail because 
they do not take into account this inhomogeneity Q . 

Modularity is one of the most important aspect of these 
networks with direct implications to transport properties. 
Here, we quantify the modular character of complex net- 
works according to our box-covering algorithm and reveal 
a connection between modularity and flow. Our results 
are consistent with recent experiments and metabolic flux 
studies, and provide a theoretical framework to analyze 
transport in a wide variety of network systems. 



A. Metabolism modeling 

In metabolism modeling, there exist three main ap- 
proaches (i) The most detailed analysis includes 
dynamic mechanism-based models Q, but in general it 
is very difficult to incorporate experimental values for 
the needed kinetic parameters, (ii) In the second ap- 
proach one simplifles the above models, and calculates 
the fluxes in a metabolic network via flux balance anal- 
ysis, which includes a family of static constraint-based 
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models The limiting factor in this analysis is that 
the problem is underconstrained (the number of unknown 
parameters, i.e. the fluxes, is larger than the number of 
metabolite conservation equations) and cannot be solved 
uniquely, (iii) Finally, a third approach that is widely 
used in metabolism modeling is to ignore stoichiometry, 
and focus only on the metabolites interactions without 
any thermodynamic aspects, which leads to interaction- 
based models jlO], i.e. undirected networks where a link 
connects two nodes that participate in a metabolic reac- 
tion. In this paper we follow this third approach and 
we use this interaction-based network to study trans- 
port on this network, by drawing an analogy between the 
metabolic network and a resistance network. If we rep- 
resent the metabolites as nodes that are linked through 
electrical resistances and the current flow represents the 
flux we can solve this problem without additional con- 
straints, and this solution may shed additional light on 
the involved processes. The advantage of our approach is 
that it can isolate the topological effect and we can ad- 
dress a broader aspect of transport in biological networks, 
such as whether the observed flux inhomogeneity is be- 
cause of the network topology itself or due to the adopted 
flux constraints. Moreover, this approach enables us to 
carry similar studies for diffusion on such networks. 



II. MODULARITY, DIFFUSION AND 
RESISTANCE 



In our work, we focus on two different examples of 
biological networks, namely the E.coli metabolic network 
and the yeast protein interaction network (PIN) [ll[ . 
Both networks have been shown to have fractal properties 
and can be covered with NsiiB) non-overlapping boxes, 
where in each box the maximum distance between any 
two nodes is less than £b, the maximum distance in a 
box [5, 6]. For a fractal network of N nodes, Nb follows 
a power-law dependence on 



NB{iB)/N ^e^'"' 
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and defines ds as the fractal (or box) dimension of a 
network. These networks are also self-similar, i.e. their 
main properties, such as the degree distribution, remain 
invariant under a renormalization scheme where each box 
is replaced by a (super) node and links between boxes are 
transferred to the nodes of the renormalized network (see 
e.g. the example in Fig. [T^ for a network G, tiled with 
= 3 boxes, that yields the network G'). Many bio- 
logical networks in the intermediate renormalized stages 
were shown to have similar properties as the original net- 
work. 

This renormalization procedure also implies the pres- 
ence of self-similar modularity in all length-scales, which 
is a central feature of these networks. The term modular- 
ity refers to the existence of sets of nodes whose links are 
connected preferably within this set rather than to the 
rest of the network. Thus, after tiling a network for a 



given value of £b , we introduce a measure of modularity 
for a network as 
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where Lf^ and L°^^ represent the number of links that 
start in a given box i and end either within or outside 
i, respectively. Large values of M correspond, thus, to a 
higher degree of modularity. Since the numerical value of 
M{£b) varies, though, with £b, a more reliable measure 
is the modularity fractal exponent dn which we define 
through: 



M{£b) ~ £i" 



(3) 



The value of dm — 1 represents the borderline case that 
separates modular (c?m > 1) from random non-modular 
{dM < 1) networks. For a lattice structure, the value of 
dM is exactly equal to dm = 1- 

In Ref. we had introduced a fractal network model 
where a network grows by adding m new offspring nodes 
to each existing network node, resulting in well-defined 
modules. In that version of the model, modules are con- 
nected to each other through x = 1 links, which leads to a 
tree structure. A generalization of this model (presented 
in detail in the Supporting Information) allows us to tune 
the degree of modularity in the network by assigning a 
larger number of links x > 1 between modules. While for 
a; = 1 all the modules are well-defined, increasing x leads 
to the presence of loops and to a progressive merging of 
modules, so that for large x values a node cannot be as- 
signed unambiguously in a given module and modularity 
is destroyed. A straightforward analytical calculation in 
this case leads to (see Supporting Information) 
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In this paper, we use this value of dj^i for the model 
and calculate dM for real networks in order to study the 
influence of modularity on network transport. 

In general, the problem of transport is expressed in 
terms of T{£; ki, ^2) and R{£; fci, /c2), where T is the av- 
erage first passage time needed by a random walker to 
cover the distance £ between two nodes with degrees ki 
and fc2, respectively, and R represents the resistance be- 
tween these two nodes. For homogeneous systems (with 
very narrow degree distribution P{k)) such as lattices 
and regular fractals, there is no dependence on ki and k2 
and the average is only over the distance £. One of the 
goals of this paper is to find the scaling of T and R in 
heterogeneous networks with a broad degree distribution 
and self-similar properties. 

In the general case of a renormalizable network, T and 
R scale in the renormalized network G" (the primes al- 
ways denote a quantity for the renormalized network) as 



T'/T = £g'^'" , R'/R ■ 



£~^ 



(5) 



The exponents dw and ^ are the random walk exponent 
and the resistance exponent, respectively. This equation 
is valid as an average of T and R over the entire system, 
applying for example in different generations when grow- 
ing or renormalizing a fractal object. Thus, this relation 
holds true for both homogeneous and inhomogeneous sys- 
tems. 

For homogeneous systems, the above exponents dw and 
C are related through the Einstein relation Q 



(6) 



where ds is the fractal dimension of the substrate on 
which diffusion takes place. This relation is a result of 
the fluctuation-dissipation theorem relating spontaneous 
fluctuations (diffusion) with transport (resistivity) and 
the underlying structure (dimensionality) Although 
the validity of this relation for scale-free networks is not 
yet clear, our following analysis shows that it also applies 
for these systems, as well. 



III. RENORMALIZATION AND SCALING 
THEORY 

The renormalization procedure on self-similar biolog- 
ical networks provides a yet unexplored method for es- 
timation of the dynamical exponents in these systems. 
Since such a network is left invariant after substitut- 
ing all nodes in a box with a single node, we can cal- 
culate the transport properties on the networks during 
successive renormalization stages. With this method we 
can also study transport in biological networks before {R 
and T) and after renormalization {R' and T'). The re- 
sults are presented in Fig. [T] for the yeast PIN, the E.coli 
metabolic network and the fractal network model with 
dM = 1-46 (tree, highly modular) and djv/ = 1-26 (net- 
work with loops and lower modularity). Figures[l]Jb) and 
[Ijc) suggest a linear relation between R' and R (T' and 
T) for a given value of £b, so that the ratio R' /R (and 
T' /T) is almost constant for all boxes in the system for 
this Ib value. For a different value of the box diameter £b 
this ratio is again constant for all boxes in the network, 
but assumes a different value. We can plot the values of 
this ratio as a function of the network size ratio N' /N 
for different values of £b (Figs. [TJi and [TJj). The data 
indicate the existence of a power-law dependence, and a 
comparison with the model networks shows that the re- 
sults are consistent with PIN exhibiting a more modular 
structure compared to the metabolic network. 

Although in principle we can use our numerical results 
to directly calculate the exponents d^ and ( through 
Eq. [5l this method is not practical because the varia- 
tion of £b is very small. We can overcome this difficulty 
by using the system size N instead, where we combine 
Eqs. (P) and d!]) to get 
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FIG. 1; (a) Example of a network G that undergoes renor- 
malization to a network G' . For this example the maximum 
possible distance between any two nodes in the same box has 
to be less than £b = 3. We compare the resistance R{A, B) 
and diffusion time T{A, B) between the two local-hub nodes 
A and B in G to the resistance R'{A',B') and T'{A',B') of 
the renormalized network G'. The values of R' and T' are 
found to always be proportional to R and T, respectively, for 
all pairs of nodes. (b),(c) Typical behavior of the probability 
distributions for the resistance R vs R' and the diffusion time 
T vs r', respectively, for a given Ib value. Similar plots for 
other Ib values verify that the ratios of these quantities dur- 
ing a renormalization stage are roughly constant for all pairs 
of nodes in a given biological network. (d),(e) We present the 
average value of this ratio for the resistance R/R' and the 
diffusion time T /T' , respectively, as measured for different 
Ib values (each point corresponds to a different value of Ib)- 
Results are presented for both biological networks, and two 
fractal network models with different (Im values. The slopes 
of the curves correspond to the exponents C,/dB in (d) and 
dw/dB in (e) (see Table 1). 
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TABLE I: Values of the exponents calculated from Fig. [T] 
Network ds C/ds d^/dn 

Metabolic network (E.coli) 3.3 0.08 ± 0.1 0.98 ± 0.1 
PIN (yeast) 2.2 0.3 ± 0.04 1.3 ± 0.04 

Model (dA/ = 1.46, m/x = 2/1) 1.46 jJ^I i -l inS 



Model (dM = 1.26, m/x = 3/2) 1.89 i ( 



In J 

r ing 

kin 2 



l(lfl + 2) 



Thus, the slopes in Figs. [TJi and [T^ correspond to the 
exponent ratios d^i/dB and C/^b, respectively. 

Notice also that the verification of the above equation 
through Fig. [T] validates the relation in Eq. ([5]) for inho- 
mogeneous systems. The numerical values for the calcu- 
lated exponents are shown in Table HI These ratios are 
consistent in all cases, within statistical error, with the 
Einstein relation, Eq. 

Using these scaling arguments and the renormaliza- 
tion property of these networks we next predict the de- 
pendence of both R and T on the distance i between 
two nodes and their corresponding degrees ki and fc2. 
After renormalization the network becomes smaller, so 
that both the degrees and the distances in the network 
decrease. A distance in G is scaled by a factor in 
G' so that £' — £/£b, while in earlier work [5] it has 
been shown that the degree k of the largest hub in a box 
transforms to a degree k' for the renormalized box, where 
k' — i'lf'^k, and dk is an exponent describing the scal- 
ing of the degree. According to the result of Fig. [1] and 
Eq. (O it follows. 
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Using dimensional analysis (see Supporting Information) 
we can show that 



R{t, ki,k2 



T{i;k,,k2)=k^-^^\fT 
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(11) 



where and /t() are undetermined functions. In 

the case of homogeneous networks where there is almost 
no /c-dependence, these functions reduce to the forms 
fnix,!) — x'' , fxix,!) — x'^^' , leading to the classical 
relations i? - ^'^ and T - £'^^ . 

The scaling in Eqs. PH)) and (fTTj) is supported by the 
numerical data collapse shown in Fig. [2l For the data 
collapse we used the values of the exponents C and dw as 
obtained from the renormalization method above (Table 
m confirming the scaling in Eqs. pUj) and pT|) . 

The functions fa and /t introduced in Eqs. pO]) and 
PTjl have two arguments, so we first need to fix the ratio 
ki/k2 and in the plot (Fig. [2]) we present different ratios 
using different symbols. We observe that the differences 
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FIG. 2: Rescaling of (a) the resistance and (b) the diffusion 
time according to Eqs. (|10|) and ([11} for the protein inter- 
action network of yeast (upper symbols) and the fractal net- 
work generation model (lower filled symbols). The data for 
PIN have been vertically shifted upwards by one decade for 
clarity. Different symbols correspond to difi'erent ratios k\/k2 
and different colors denote a different value for fci. Inset: Re- 
sistance _R as a function of distance i, before rescaling, for 
constant ratio k\/k2 = 1 and different k\ values. 



among varying ratios are small, so that the fci/fc2 depen- 
dence in Eqs. (fTO|) and pT|) can be neglected. 



(12) 



k 



with an analogous form for T{£; ki, k2)- We are, thus, 
led to a simpler approximate form where the values of 
diffusion and resistance between two nodes depend only 
on the lowest degree node ^2 and the nodes distance £. 
This is a generalization to real networks of the result 
in Ref. |12|| for random model networks where, without 
taking into account the distance, the resistance was found 
to depend solely on k2- 



IV. INFLUENCE OF MODULARITY ON 
TRANSPORT 

Modularity is a central feature of biological networks 
which contributes to a more efficient use of resources in 
the network, with yet unclear consequences for transport 
in these systems. In this section we use the fractal net- 
work model, which reproduces the main features of real 
networks, to better understand the influence of modular- 
ity on transport. 

A direct calculation of ( for the fractal network model 
is as follows. We consider the growth model where the 
distance £ between two nodes increases by a factor of 3, 
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i.e. I' = 3€, and the resistance R between two neighbor 
nodes in G increases by a factor 3/x in G", since the 
hnear distance between the two nodes has increased by 
a factor of 3 and there are x parallel paths connecting 
these nodes, i.e. R' — 3R/x. Combining these equations 
with Eq. ^ wc find that the exponent C, for this model 
is given by: 

C = i!^ = i-i!i|. (13) 



In 3 



In 3 



1 



Notice that for a tree structure (x = 1) we get ( 
as expected. This result is also an important step in 
linking statics with dynamics (a long standing problem 
in percolation theory [1, Using only the value of x 
that describes how self-similar modules are connected to 
each other we can directly obtain the dynamic exponent 
C, i.e. how the structural property of modularity affects 
dynamics. 

For the fractal dimension of model networks we already 
know that M 



dn = 



ln(2m + x) 

hTs 



(14) 



If we assume that the Einstein relation in Eq. ([6]) is valid, 
then we can also calculate the value for the random walk 
exponent dw' 

d.-^tr^- (15) 

A comparison of Eq. with Eq. (|4|) yields: 

dw = 1 + du ■ (16) 

This is a very simple, yet powerful result. It manifests 
that for the fractal network model the degree of mod- 
ularity directly affects the efficiency of transport and is 
the main feature that controls the type of diffusion. 

The above relation, Eq. p^ . is verified in Fig. [3] for 
the fractal network model. We generate a number of 
model networks where we vary both the number of loop- 
forming links x and the number of offsprings m. For each 
pair of m and x we calculate numerically the exponent 
dm from the slope of figures similar to Figs. [TJi and [TJ; 
and use Eq. ^ for the value of dM- The results are 
fully consistent with Eq. pB|) and all the points lie on 
the predicted line. Sub-diffusion {d^ > 2) is observed 
for dM > 1, in accordance with our observation that 
modularity slows down diffusion. On the contrary, for 
non-modular networks diffusion is accelerated remark- 
ably (c?tt, < 2) which is also in agreeement with previous 
work on random networks. When d^j — 1 we recover 
classical diffusion (dw — 2), even though the structure is 
still that of a scale-free network. 



V. FLOW DISTRIBUTION ACROSS THE 
NETWORK 

In our scaling theory above we derived results for the 
average values of the current flowing in a complex net- 
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FIG. 3: Comparison of the random walk exponent ex- 
tracted numerically (symbols) with the theoretical prediction 
(Eq. (|16|l . line) vs the modularity exponent dM, for different 
values of m and x. Inset: Direct (unsealed) numerical calcu- 
lation of dw as a function of m, for varying x values (shown 
in the plot). 



work. The inherent inhomogeneity and modularity of 
biological networks is expected, though, to strongly influ- 
ence the distribution of flow throughout the network. Us- 
ing flux balance analysis, it was recently shown that the 
distribution of fluxes in the metabolic network is highly 
uneven and a small number of reactions have the largest 
contribution to the overall metabolic flux activity [l| . To 
study the influence of the complex substrate on the flow 
distribution we calculate the probability P{I) of current 
lij flowing between all two neighboring nodes i and j . 

The probability distribution P{I) for the magnitude 
of the current / across a link of the metabolic network 
decays according to a power law (see Fig. 2]). The form 
of the curve and the exponent in the range 1.0-1.5 of the 
decay is very similar to those found in previous studies 
of the metabolic flux, both experimental fl^ and theo- 
retical 1] . The decay suggests that only a small fraction 
of link carries high current. For the yeast PIN the distri- 
bution is even broader, and its form is different from the 
metabolic network. The variation of P{I) is smaller in 
PIN compared to the metabolic network, meaning that 
in PIN there is a larger number of important links that 
carry large currents. The self-similar character of the bio- 
logical networks is also verifled for the distribution P(/), 
as well, which remains invariant after renormalizing the 
network, as seen in the inset of Fig. U) 

Next we explore the connection between P(/) and 
topology. The information contained in P{I) can be bet- 
ter understood if we compare these results with a sur- 
rogate random case. The random case is obtained by 
rewiring the original networks, preserving the degree dis- 
tribution P{k), but destroying any correlations between 
neighboring nodes. Thus, we remove all traces of the ini- 
tial network organization. The distribution P{I) for the 
metabolic network remains almost the same under this 
rewiring, indicating that, despite the modular character 
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FIG. 4; Probability distribution P{I) of current magnitudes 
I flowing through the links in PIN (solid triangles) and 
metabolic networks (solid circles). Empty symbols are the 
corresponding results for the randomly rewired networks. In- 
set: Invariance of P{I) for the metabolic network under renor- 
malization with different is values. 



of the network, the original structure behaves similar to 
uncorrelated networks, since the degree correlations do 
not affect P{I) in the metabolic network. In contrast, 
the distribution P{I) for the rewired PIN is very differ- 
ent than the original distribution and is similar to that of 
the metabolic network. We can, thus, conclude that the 
original PIN has a much richer structure that deviates 
from the random case, corroborating thus the results on 
modularity from the previous section on diffusion. 

The above results for the biological networks can be 
understood in terms of the fractal network model, where 
we can control the number of links that form loop struc- 
tures in a network. In Fig. [5] we calculate P{I) for the 
model network with to = 4 and x = 2, which is a highly 
modular structure {dM = 1.46). The form of the P{I) 
distribution is similar to that of the PIN, but if we add 
a small number of random links (or equivalently rewire a 
small part of the network) this distribution is significantly 
influenced in a similar way as observed in random PIN 
rewiring. This suggests again that in the case of PIN 
modularity is high. In the inset of this Figure we can 
also see that as x increases, i.e. more loops appear in the 
structure, the distribution has a longer tail, which shows 
that there is a smaller number of high-current links. 

Since the number of added links in Fig. O is small, the 
modularity is preserved. We verified that the main fac- 
tor that influences P{I) is the number of loops in the 
network, rather than modularity itself, by fixing dM and 
only varying x. In this case (described in Supporting 
Information) the P{I) distribution is different for net- 
works with the same d^f exponent. This can also be seen 
through Eq. (fT^]) , where the resistance exponent depends 
only on x. 

Using the information of we can also construct the 
'backbone' of the network in the form of the minimum 



FIG. 5: Probability distribution P{I) for the fractal model 
before and after randomly adding 1% of links or rewiring 10% 
of the network. Inset: P{I) for the fractal model with varying 
dM values, where m — 2 and x varies from 2 to 4. 

spanning tree (MST) [3]. The importance of such a 
tree is that it identifies the substructure of the network 
that is dominant for transport. Starting from a com- 
pletely empty network we insert links in decreasing order 
of current magnitude, provided that they do not form a 
loop. The resulting MST tree for the PIN is presented 
in Fig. [6l where the thickness of the links in the drawing 
increases logarithmically with increasing current through 
a link. The color of a node corresponds to the function 
performed by a protein in the network. This tree (created 
solely on the base of current flowing through a network) 
exhibits a large degree of modularity where nodes that 
perform similar functions are close to each other. It is 
also possible through this construction to identify the 
most critical links in the network in terms of the largest 
current flowing through them. 

Thus, the emerging picture from the above analysis 
for the PIN is one of a network with a strong backbone 
that carries most part of the flow combined with loops 
organized mainly within modules, so that flow through 
this backbone is not really influenced. This result high- 
lights the strong modularity in the PIN structure. If a 
structure has a smaller degree of internal organization, 
as is the case in the metabolic network, then flow is more 
uniformly distributed. 



VI. SUMMARY 

Summarizing, we have presented scaling arguments 
and simulations on a class of self-similar complex net- 
works, concerning transport on these networks. Diffusion 
and resistance in these biological networks is important, 
since they are both estimates of how many paths con- 
nect two nodes and how long these paths are. By using 
the self-similarity property and a network renormaliza- 
tion scheme we have developed a scaling theory of the 
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FIG. 6: (a) Current flow through the Unks of the yeast PIN 
network, for one random selection of the two nodes acting 
as current input/output, (b) Minimum spanning tree for the 
PIN. The thickness of a link corresponds to the current flow- 
ing through this link. Different node colors correspond to 
different protein functions. 

resistance and diflFusion dependence on both the distance 
between two nodes and their corresponding degree. We 
were able to recover a relation between network modu- 
larity and transport, while the flow distribution in these 
networks was found to be consistent with earlier studies, 
using different approaches. 

VII. METHODS 

A. Resistance measurements 

In order to measure the conductivity between two 
nodes A and B we consider that all links in the underly- 
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ing network between any two neighbor nodes i and j have 
unit resistances Rij = 1. By fixing the input current to 
//I = — 1 and the output to = 1 we can solve the 
Kirchhoff equations and compute the voltages in the sys- 
tem. The measured resistance is then Ra-^b — Va — Vb- 
However, due to the required inversion of the relevant 
matrices we are limited by the computer resources to 
networks of relatively small size, i.e. N < 10'' nodes. 

In principle, the magnitude of /y depends on the se- 
lection of the current input /output nodes. Upon closer 
inspection, though, we found that the distribution of the 
current magnitudes in the network links is not very sen- 
sitive to the selection of the current source and the cur- 
rent sink. Comparison of the result of averaging over one 
input and many outputs and over twenty different in- 
put and output pairs of nodes for the metabolic network 
showed that within statistical error these two results are 
almost identical. 



B. Diffusion measurements 

In many cases (and especially those including real- 
life networks) direct measurements of diffusion on com- 
plex networks exhibiting the small-world property may 
present significant difficulties, due to the limited time- 
range where diffusion takes place before settling quickly 
to a distance equal to the typical (very short) network 
diameter. The rising part of the mean-squared displace- 
ment as a function of the time is very small and reliable 
measurement of the diffusion exponent is very hard to 
do. Moreover, we need additional information in order 
to quantify the k and i dependence of the diffusion time 
between two nodes. For this purpose we used the peak of 
the first passage time distribution as a typical diffusion 
time T{A, B) between two points A and B in the net- 
work. Since this quantity may be assymetric depending 
on which node we consider as origin, our diffusion time 
T represents the average of {Ta^b + Tb^a)/"^- 
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